Nucleic Acids Research

Transfer RNA (tRNA) molecules serve as essential adapters during protein translation. While direct RNA sequencing (DRS) via Oxford Nanopore Technologies has emerged as a powerful platform for systematic tRNAome profiling, we currently lack a simple and robust statistical framework for nanopore tRNA data analyses. Here, we address this gap by developing AMaNITA (Abundance, Modifications, and Nanopore Intensity Toolbox Application), an end-to-end bioinformatic workflow that enables simplified, robust, and scalable analyses of nanopore native tRNA sequencing datasets. AMaNITA streamlines the entire analytical trajectory: from upstream processing (basecalling, mapping, filtering, batch effect correction) to downstream assessment of differential tRNA abundance and modification stoichiometry. The workflow generates an interactive HTML report for data exploration and analysis, allowing the user to download the source data files and resulting plots. AMaNITA can be executed using Singularity from the command line, without requiring installation of dependencies.

Deinococcus radiodurans harbors a largely classical bacterial DNA repair machinery yet displays exceptional resistance to ultra-violet and ionizing radiation. To investigate whether crosstalk between its DNA repair pathways contributes to this phenotype, we mapped putative interactions between the nucleotide excision repair (NER) and base excision repair (BER) pathways, which together are responsible for the removal of nucleobase lesions. Using a bacterial two-hybrid system, we identified multiple direct interactions between NER and BER proteins, notably involving the two UvrA variants, and validated these interactions in vitro. Furthermore, functional analyses revealed that NER interferes with the BER-mediated removal of oxidized guanines by the Fpg DNA glycosylase, likely through competition for DNA binding and sequestration of Fpg. Finally, UvrB and UvrC were found to further process the Fpg incision product in an ATP-dependent, UvrA1-independent manner. Together, these results demonstrate a multi-level crosstalk between NER and BER in D. radiodurans, which may contribute to its extraordinary DNA repair capacity. To our knowledge, this represents the first evidence of such a complex interplay in bacteria.

Many bacteriophages encode anti-CRISPR (Acr) proteins that inhibit the CRISPR-Cas immune systems. Rapid acr gene expression upon phage entry enables CRISPR-Cas neutralisation, but can impact phage fitness if unregulated. Therefore, Acr production is often controlled by distinct families of co-encoded anti-CRISPR-associated (Aca) proteins, which are usually helix-turn-helix (HTH) regulators that bind DNA within acr-aca operon promoters. Previously, we demonstrated that the Aca2 family additionally represses Acr production translationally by binding structured RNA motifs within the 5' untranslated regions (UTRs) of the acr-aca mRNA. Here, through systematic bioinformatic analyses, we provide evidence of structured RNA motifs in the 5' UTRs of operons encoding members of other Aca families, and that Aca1 also specifically binds its cognate RNA motif. Additionally, many Aca proteins are predicted to regulate not only their own but also adjacent operons with potential anti-defence genes. Indeed, we show that Aca14, newly identified in this study, represses two predicted anti-defence operons. Aca14 is a ribbon-helix-helix domain protein, revealing regulatory diversity beyond the canonical HTH Aca family members. Collectively, our findings expand the understanding of acr regulation in mobile genetic elements and reveal novel mechanisms by which phages fine-tune anti-defence gene expression. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/740549v1_ufig1.gif" ALT="Figure 1"> View larger version (21K): org.highwire.dtl.DTLVardef@1ca8b1dorg.highwire.dtl.DTLVardef@73291eorg.highwire.dtl.DTLVardef@90ab45org.highwire.dtl.DTLVardef@1dc739e_HPS_FORMAT_FIGEXP M_FIG C_FIG

Artificially Expanded Genetic Information Systems (AEGIS) increase the information content of nucleic acids by including new nucleobase pairings that are orthogonal to those of canonical Watson-Crick nucleobases. DNA ligases do not form direct interactions with the nucleobases during catalytic turnover, suggesting that these enzymes should efficiently and faithfully join double-stranded AEGIS substrates. Here we report the systematic investigation into the validity of this hypothesis for structurally-diverse DNA ligases employing substrates built from the eight nucleotide hachimoji genetic alphabet, where orthogonality is achieved by rearranging the hydrogen bonding patterns seen in canonical Watson-Crick pairs. We find that single, or multiple, non-canonical bases are well tolerated at the 5 prime-end of the nick. However, tracts of consecutive non-canonical bases at the 3 prime-end of the break significantly decrease ligation efficiency or abolish it altogether. Possible reasons for this apparent bias against non-canonical nucleobases could include incompatibility in electrostatic interactions between the ligase active site and the non-canonical substrates or altered conformational preferences and/or dynamics in key catalytic intermediates. We also observe single hachimoji mismatches are ligated more frequently than mis paired canonical bases, potentially due to promiscuous pairing of tautomeric forms of the non-canonical bases.

Transcriptional termination efficiency is considered an important parameter for fine tuning bacterial gene expression. Still, the design principles that determine transcription termination efficiency remain poorly understood. In this study, we aimed to investigate the impact of the 3' untranslated region (3'UTR) on gene expression in Escherichia coli and other bacteria. First, 3'UTR variant sequences were generated, with randomized 30 bp sequences inserted between the STOP-codon and an intrinsic terminator, consisting of a GC-rich hairpin and a downstream poly(U)-tail. Using three reporter genes, it was found that different 3'UTR sequences resulted in an up to five-fold difference in protein production, independent of the upstream coding sequence. The highest protein production was achieved when an adenosine was present directly upstream of the terminator hairpin. This was consolidated by systematic substitution of key nucleotides of the terminator and assessing their effect on mRNA and protein levels. Subsequently, we developed a predictive random forest machine learning model trained on the termination efficiency of different natural and synthetic terminator sequences, revealing an important role for the nucleotides directly upstream of the terminator hairpin. Altogether, this study showed that an additional adenosine nucleotide upstream of the terminator hairpin leads to improved protein production while reducing terminator read-through.

The MSH1 gene in Arabidopsis thaliana (AtMSH1) encodes a modular enzyme that consists of an N-terminal MutS DNA mismatch repair module fused to a C-terminal GIY-YIG nuclease. Disruption of MSH1 reverses the low organellar mutation rates that distinguish plants from other eukaryotes. However, the precise mechanism by which MSH1 prevents the accumulation of mutations remains unclear. Here, we show that AtMsh1 accurately recognizes and cleaves dsDNA containing mismatches and short indels. AtMsh1 efficiently cleaves dsDNA containing lesions generated by oxidative damage or deamination, with a strong preference for U:G mismatches. AtMsh1 cleaves DNA through an ATP-dependent enzymatic mechanism that requires divalent metal cofactors such as Mg2+. The enzyme introduces incisions at defined positions relative to the lesion or mismatch: approximately nine nucleotides 5' of the mismatch on the affected strand and twelve nucleotides 3' on the complementary strand. This offset cleavage generates staggered DNA ends with three-nucleotide overhangs. Although AtMsh1 displays positional specificity in its cleavage activity on substrates containing lesions and mismatches, it exhibits nonspecific double-stranded DNA cleavage in the presence of Mn2+. These findings establish AtMsh1 as a minimal mismatch repair (MMR) system in which mismatch/lesion recognition and DNA cleavage are functionally coupled. We propose that the resulting dsDNA breaks are processed by exonucleases that mediate single-stranded DNA resection, thereby removing the mismatch or lesion while generating a 3' single-stranded DNA overhang suitable for homologous recombination (HR) repair and gene conversion. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=82 SRC="FIGDIR/small/731605v1_ufig1.gif" ALT="Figure 1"> View larger version (17K): org.highwire.dtl.DTLVardef@156ce6eorg.highwire.dtl.DTLVardef@abb52org.highwire.dtl.DTLVardef@90286aorg.highwire.dtl.DTLVardef@4d7dce_HPS_FORMAT_FIGEXP M_FIG C_FIG Mismatch recognition and nuclease activity by plant organellar MutS Homolog 1 drive organellar genome maintenance.

Site-directed RNA editing (SDRE) utilizing the adenosine deaminase acting on RNA (ADAR) enzymes is commonly facilitated by guide RNAs (gRNAs) optimized to enhance on-target editing and minimize bystander effects. However, the impact of gRNA binding and ADAR-mediated SDRE on canonical pre-mRNA splicing remains poorly understood. Here, we developed an in vitro transcript-specific editing strategy that enables selective targeting and direct comparison of SDRE in pre-mRNA and mature mRNA. Using splice-relevant variants associated with inherited retinal diseases, we investigated the effects of SDRE on exonic, near-canonical intronic, and deep intronic splice variants. We identified gRNA-induced splice perturbation at exonic and intronic targets and observed that higher editing levels could be associated with increased splice disruption. Conversely, SDRE of two exonic splice variants and a deep intronic variant resulted in increased production of correctly spliced transcripts, demonstrating the potential of SDRE for splice modulation. Finally, by dissecting the effects of ADAR expression and gRNA design on editing and splicing outcomes, we established a system for identifying design principles that reduce splice interference and enhance the generation of correctly spliced, edited transcripts. These findings highlight the importance of considering transcript context and splicing consequences in the development of SDRE-based therapeutic strategies. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=92 SRC="FIGDIR/small/741163v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@2e4703org.highwire.dtl.DTLVardef@3834d3org.highwire.dtl.DTLVardef@1409c8org.highwire.dtl.DTLVardef@18e463b_HPS_FORMAT_FIGEXP M_FIG Created in BioRender. Schneider, N. (2026) https://BioRender.com/59i0tiy C_FIG

During replication, nucleolytic processing of apurinic/apyrimidinic (AP) sites in single-stranded (ss)DNA is attenuated by the evolutionarily conserved 5-hydroxymethylcytosine binding, embryonic-specific (HMCES) protein. HMCES forms a covalent thiazolidine linkage with the ring-opened aldehyde form of a ssDNA AP site to stabilize the AP site and suppress the formation of DNA double-strand breaks. The resulting HMCES DNA-protein cross-link (DPC) can then be digested by the SPRTN protease and bypassed through mutagenic translesion synthesis (TLS). Here, we use Xenopus egg extracts and molecular dynamics simulations to investigate how HMCES-DPC formation influences the mutagenicity of AP site bypass. We show that SPRTN processes the HMCES-DPC to a five amino acid peptide adduct prior to TLS. Surprisingly, the mutagenicity of HMCES-DPC bypass is insensitive to the extent of DPC proteolysis and depends only on cross-link formation, which traps the AP site in a more dynamic ring-opened configuration. We further show that the spectrum of mutations produced during bypass of HMCES-adducts strongly depends on the template strand nucleotide immediately 5 of the AP site. Our data support a model in which HMCES-DPC formation increases the conformational flexibility of the DNA template, allowing template slippage and use of the 5 template nucleotide to direct insertion opposite the adducted AP site.

Various methods have detected miRNA-target interactions via immunoprecipitation of UV-crosslinked Argonaute ribonucleoprotein complexes, followed by intermolecular ligation of bound miRNAs to target strands, forming chimeric RNAs. To date, these methods have relied on conventional viral reverse transcriptases (RTs) to generate cDNAs for sequencing. However, crosslinked RNAs often retain adducts after purification, which can make them poor templates for viral RTs. Here, we adapted OTTR (Ordered Two-Template Relay) techniques to generate cDNAs from Ago2-bound RNAs. OTTR makes use of a modified retroelement-encoded RT, which is strongly processive even on templates with modifications or adducts. We show that this "OTTR-CLASH" method increases the frequency of generating chimeric RNAs compared to previous methods. We also developed an improved bioinformatic pipeline for analysis of these data, and we use this to catalog miRNA-target interactions not previously described in the literature. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=147 HEIGHT=200 SRC="FIGDIR/small/738487v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@13bc276org.highwire.dtl.DTLVardef@5beb41org.highwire.dtl.DTLVardef@b204e5org.highwire.dtl.DTLVardef@15f747d_HPS_FORMAT_FIGEXP M_FIG C_FIG

Abstract The universe of possible nucleotide sequences expands combinatorially with sequence length, vastly exceeding the fraction sampled by real genomes. Yet genomic sequences exhibit reproducible compositional symmetries and recurrent structural motifs, indicating that biological sequence space is shaped by strong organizing constraints. Here, we introduce an explicit framework for constructing and visualizing the complete sequence universe using the Newtonian polynomial for a four-letter alphabet, and for identifying biologically relevant subsets through the application of fundamental filters. Three filters of biological relevance are formulated: (i) the constraint that DNA predominantly exists as an antiparallel-stranded double helix, (ii) the second Chargaff parity rule, which enforces approximate strand symmetry in single-stranded sequence composition, and (iii) genome shadows, reflecting the imprint of concerted sequence changes. Successive application of these filters dramatically reduces the accessible sequence space and reveals distinct symmetry classes. Among these, mirror-symmetric sequences occupy a privileged position because they are invariant under strand reversal and therefore compatible with both antiparallel and parallel strand orientations. This dual compatibility enables such sequences to bridge otherwise disjoint structural subspaces of DNA. G-rich members of this class are shown to have a strong propensity to form G-quadruplex architectures that incorporate parallel-stranded domains while remaining compatible with duplex DNA. We propose that this structural versatility provides a mechanistic basis for the recurrent association of G-rich mirror-symmetric sequences with recombination hotspots and genome rearrangements. Together, these results establish a symmetry-based framework for understanding how combinatorial sequence space is filtered into biologically functional DNA motifs.

In the evolving landscape of RNA research, the classification and analysis of RNA motifs is necessary to uncover the intricate mechanisms governing cellular and viral processes. Here we apply the coarse-grained RAG (RNA-As-Graphs) framework to advance the classification and understanding of RNA motifs, with a focus on expanding the RNA Motif Atlas through the inclusion of novel viral RNA structures. By analyzing 273 experimentally resolved viral RNA structures from the Protein Data Bank (PDB) using RAG dual-graph representations, we identify 14 previously uncatalogued viral RNA motifs. These motifs, which include tRNA mimicking domains, exoribonuclease-resistant domains, and internal ribosome entry sites, expand the diversity of RNA to a total number of 197 dual graph motifs. We applied k-means, PAM, and Ward clustering and observed substantial overlap between viral and general RNAs. The expanded library of RNA motifs and submotifs provides a resource for motif discovery and RNA design.

Analysis of RNA-binding proteome data from different bacterial species revealed many cell wall metabolic enzymes cross-linking to RNA in vivo, hinting that these proteins directly bind RNA. Surprisingly, penicillin-binding proteins (PBPs) were also abundantly identified as putative RNA-binding proteins. The cell surface localisation properties of many of these proteins therefore beg the question at what stage of their cellular life cycle these proteins interact with RNA and what the functional significance is. Here, we characterised the RNA-binding activity of PBP2a, the alternative transpeptidase that confers {beta}-lactam resistance in MRSA. Using in vivo RNA-binding assays, we show that PBP2a interacts with hundreds of transcripts without apparent sequence specificity. Computational analyses identified a possible RNA-binding cleft in PBP2a proximal to its active site. Mutation of only two predicted positively charged residues located in this cleft substantially reduced cross-linking in vivo, implying that RNA recognition is largely dictated by RNA backbone interactions. While PBP2a does not regulate RNA steady-state levels, RNA-binding appears important for proper protein function: an RNA-binding deficient mutant exhibits reduced oxacillin resistance. These findings establish PBP2a as an RNA-binding protein in vivo and provide a framework to investigate how this non-canonical interaction may relate to cell wall biogenesis and {beta}-lactam resistance.

RNA-based medicines rely on modified nucleotides to promote immune evasion and in vivo efficacy. Nucleotides generated from RNA degradation are either exported or recycled through metabolically favorable salvage pathways, though whether modified nucleotides are efficiently recycled remains unclear. N4-acetylcytidine (acC) is a naturally occurring modification in rRNA and tRNA that has shown promise in therapeutic mRNA applications. However, N4-acetylation impairs cytidine deamination, the first step in cytidine salvage. Here, we investigate the endogenous mechanisms that enable acC metabolism. Through sensitive sequence and structural analyses, we identify the uncharacterized human ASCH domain protein EOLA1 as a key acC deacetylase in nucleotide salvage. EOLA1 inactivation leads to free intracellular acC accumulation and increased cytotoxicity upon nucleotide export inhibition. While steady-state acC levels in cellular RNAs remain unchanged, EOLA1-dependent regulation of free acC is evident basally and is exacerbated by exogenous mRNA delivery. Proteomic analyses place EOLA1 in proximity to ribosomal proteins, adjacent to endogenous acC sources. In vitro assays confirm EOLA1 specificity for acC, and structural analysis reveals a narrow nucleotide-binding pocket consistent with mononucleotide selectivity. These findings identify EOLA1 as a bona fide acC eraser and uncover a previously unrecognized pathway for recycling modified nucleotides with relevance to therapeutic RNA design.

The ribosome, long regarded as a passive, uniform machine, has only recently been recognised as a direct regulator of translation. Mass spectrometry and sequencing approaches have shown that heterogeneity in ribosome composition exists, which can actively regulate the translational process. One source of this heterogeneity is the modification of ribosomal RNA (rRNA), primarily pseudouridylation (pseU) and 2'-O-methylation (2OMe), mediated by specific H/ACA and C/D box small nucleolar RNAs (snoRNAs). Here, we investigate how the stoichiometry of rRNA modifications varies during embryonic stem cell differentiation. Using the modification basecalling capability of Nanopore direct RNA sequencing, we have identified distinct stoichiometric changes in modification patterns between pluripotent and differentiated cells, revealing highly dynamic, site-specific regulation. Further, profiling of snoRNA expression during trilineage differentiation revealed differential expression of H/ACA and C/D box snoRNAs responsible for a subset of these dynamic modifications. By integrating rRNA and snoRNA sequencing approaches, we have built a comprehensive profile of rRNA modification dynamics during early embryonic cell fate decisions, highlighting potential regulatory mechanisms for ribosome heterogeneity during development. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=191 HEIGHT=200 SRC="FIGDIR/small/743918v1_ufig1.gif" ALT="Figure 1"> View larger version (26K): org.highwire.dtl.DTLVardef@48ca91org.highwire.dtl.DTLVardef@eb0418org.highwire.dtl.DTLVardef@159fc8corg.highwire.dtl.DTLVardef@d34f19_HPS_FORMAT_FIGEXP M_FIG C_FIG

MYC oncoprotein binding at promoters and enhancers influences RNA polymerase II (RNAPII)-driven gene expression. Numerous genes also bind MYC near their transcriptional end sites (TESs). This often allows direct promoter-TES contact via looping and further regulates total and 'read-through' transcription that extends beyond standard termination sites. We aimed here to better clarify the rules governing TES associated MYC and/or RNAPII binding cross-talk in human and murine cells. Using ChIPseq and RNAseq datasets from the ENCODE portal and elsewhere, MYC and RNAPII binding profiles were found to differ around TESs and transcriptional start sites (TSSs). Variations in E box flanking sequences likely accounted for the somewhat lower affinities of MYC for TES-associated sites. Motifs for numerous other transcription factors were also observed to cluster non-randomly and in close proximity to MYC and RNAPII binding site peak summits. On average, genes with TES-proximal MYC or RNAPII sites were more highly expressed than those without, although co-binding tended to be suppressive. Both normal and neoplastic proliferative stimuli altered the MYC and RNAPII binding patterns of many genes, indicating that 'category switching' was common, subject to disparate external signals and often reversible. Functionally related gene sets with high levels of read-through transcription were uniformly marked by significant amounts of TES-associated MYC and/or RNAPII binding. These findings indicate that, both independently and together, MYC and RNAPII binding near TESs dynamically impact total and read-through transcription while also coordinating the expression of many common purpose gene sets.

Small regulatory RNAs (sRNAs) are key players in bacterial adaptation to stress. They often occupy central positions in regulatory networks and control the expression of multiple targets. In a striking example of this, the enterobacterial OmrA and OmrB paralogous sRNAs are known to regulate about ten different targets, with extensive data suggesting the regulon is in fact much larger. Here we performed transcriptome and proteome analyses and identified more than fifteen new targets of Escherichia coli OmrA and OmrB. We validated several, including genes involved in central carbon metabolism and fatty acid synthesis, among which ppc, actP and fabA. Consistent with a role in carbon metabolism, overproducing OmrA or OmrB inhibited growth on glucose minimal medium. The analysis of suppressor mutants shows that this is due to a decreased carbon flux through the TCA cycle. Incorporating other datasets such as RIL-seq, we generated a multi-omics-based prediction of target candidates. Together, our results show that OmrA/B base-pair to various regions of their mRNA targets, and therefore likely act through diverse regulatory mechanisms. Hence, this work extends the OmrA and OmrB regulons, establishes an unsuspected connection with carbon usage, and shows the benefits of combining global analyses to investigate sRNA regulons.

DNA replication achieves error rates as low as 10-9-10-11 per base pair through the combined contribution of base selection, exonuclease proofreading, and mismatch repair. Among those processes, base selection determines the initial level of accuracy and is strongly modulated by local sequence context. Existing models address this dependence either by fitting individual rate constants for each sequence context or by invoking the global template properties, neither of which derives sequence dependence from the underlying thermodynamics and kinetics of base pair formation. Here we present a mechanism for sequence-dependent base selection fidelity, built from two physical properties: nearest-neighbor stacking thermodynamics and directional kinetic asymmetry. The model fits the experimentally observed mutation spectra from three mismatch repair-deficient organisms well (r=0.74, 0.70, and 0.63), and predicts that base-selection accuracy varies non-monotonically with temperature in a sequence-dependent manner. Our model, therefore, provides a framework that connects sequence-dependent thermodynamic and kinetic effects during nucleotide incorporation to experimentally observed mutation rates.

Translesion DNA synthesis (TLS) enables DNA replication across damaged DNA and promotes stress-induced mutagenesis that contributes to antibiotic resistance in bacteria. The conserved ImuABC mutasome is essential for TLS in many bacterial species, yet the molecular function of its accessory protein, ImuA, has remained elusive. Here we show that Mycolicibacterium smegmatis ImuA assembles into a hexameric complex, likely arranged as a dimer of trimers, with dual enzymatic activities that reshape current models of its role in DNA damage tolerance. We show that ImuA functions as an ATP-dependent helicase that preferentially unwinds DNA substrates containing single-stranded DNA overhangs and identify amino acids required for both hexamer formation and helicase activity. Unexpectedly, ImuA also possesses ATP-independent 5' exonuclease activity, selectively processing ssDNA substrates with free 5' ends. We show a basic patch on the N-terminus is essential for stabilizing both the nuclease motif and oligomerization. Together, these findings identify ImuA as an active DNA-processing enzyme rather than a passive accessory factor and establish oligomerization as a prerequisite for its function. Our work provides a mechanistic framework for understanding how ImuA may function within the ImuABC mutasome to coordinate DNA processing during translesion synthesis.

UvrD-family SF1A helicases play a variety of biological roles in DNA metabolism including replication, recombination, DNA repair, and conjugative plasmid transfer. The family is described by a subdomain architecture consisting of two RecA motor subdomains (1A and 2A) that each contain an auxiliary B-domain insertion (1B and 2B). Monomeric UvrD-family enzymes possess ATPase and 3-5 single-stranded DNA translocase activity but lack helicase activity. An enzyme dimer formed through a 2B-2B domain interface acts as a processive DNA helicase. A partial explanation for this observation, based on structural comparison between monomeric and dimeric DNA bound complexes, is that dimerization removes inhibitory contacts between the 2B domain and the double-stranded DNA, thus activating the helicase. However, biochemical observations reveal additional aspects of the dimeric mechanism that cannot solely be explained by movement of the 2B subdomain. Here, we present data showing that the Mycobacterium tuberculosis UvrD1 dimer interacts with the displaced strand of the duplex (i.e., the 5-3 strand in the direction of unwinding). In particular, bulky lesions on the displaced strand-including a thymine dimer-led to more efficient unwinding over a finite range of duplex lengths. These findings suggest a model for dimeric UvrD-family unwinding wherein one subunit makes intimate contacts with the displaced strand and that this interaction promotes the processive DNA unwinding unique to dimers.

Advances in high-throughput sequencing have associated millions of putative genetic variants with disease. However, scalable experimental methods to establish causal relationships between genetic variants and downstream transcriptional outcomes remain a major challenge. Single-cell methods that integrate genotyping with transcriptomic profiling provide a way to address this, but do not enable pre-sequencing enrichment of correctly edited cells, limiting scale. We present SELECT-seq (SNP Enrichment Leveraging Cas12a Targeting), a rapid method that allows SNP-specific PCR amplification and Cas12a-mediated fluorescence detection simultaneously with whole-transcriptome amplification. This one-pot workflow enables identification and enrichment of SNP-bearing single cells, making a rapid and scalable methodology for analysis of genotype-phenotype linkage avoiding laborious single cell cloning steps. As a proof of principle we show that SELECT-seq distinguishes U-2 OS and T-47D cell lines based on a PIK3CA (NM_006218.4:c.3463A>G) mutation while preserving transcriptome integrity. It physically enriches a rare NRF2 T80K (NM_006164.5:c.390C>A) mutant cells (6.7%) from a prime-edited pool, achieving 86% genotype accuracy, and shows 87.5% directional concordance in the transcriptomic effects compared with a clonal NRF2 T80K cell line. SELECT-seq thus provides a rapid, scalable and widely accessible approach for mapping genotype-phenotype relationships at single-cell resolution.